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I. INTRODUCTION 



Neutrino oscillation experiments have demonstrated that neutrinos are massiveparticles 
With the most recent experimental data by the MINOS and KamLAND [s] collab- 
orations atmospheric and solar mass-squared differences are now Icnown very precisely and 
global fits to all neutrino oscillation data ^ also give quite accurate determinations for the 
corresponding neutrino mixing angles. For the overall mass scale of neutrinos and the third 
neutrino mixing angle currently only upper limits exist, but considerable progress is expected 
from future double beta decay |5| and reactor neutrino oscillation experiments. 

Neutrino masses provide the first experimental signal of physics beyond the standard 
model (SM). From an experimental point of view, neutrino oscillation data can easily be 
fitted in very much the same way as the SM accounts for quark masses and mixings, i.e. 
namely by Dirac neutrino masses. From a theoretical point of view, however, such an ansatz 
is ad hoc since, being electrically neutral, neutrinos are expected to be Majorana particles 
Indeed, as noted already in j^, the dimension-5 operator, 
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{HL){HL), 



(1) 



induces Majorana masses for neutrinos once the electroweak symmetry breaks. This way 
the smallness of the neutrino masses can then be attributed to the existence of some lepton 
number violating scale larger than the electroweak scale. A variety of ways to generate this 
operator have been suggested. The resulting Majorana neutrino masses can be suppressed 
either by loop factors, by a large mass scale, by a small scale whose absence enhances the 
symmetry of theory, or by combinations of these mechanisms |lol |. 



Electroweak scale models, such as, for example the Zee model 



12] , supersymmetric models with violation of R-parity 



n 

161 



the Babu-Zee model 



171 ] or lepton number 



violating leptoquark models [18|] generate neutrino masses at loop-level, resulting in / ^ 1 
and A need not be much larger than mw- A similar situation arises in models like the inverse 
seesaw [l^. Such low-scale models have the advantage that the new fields responsible for the 
generation of neutrino masses may be directly accessible to future accelerator experiments, 
see for example 
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The most popular mechanism to generate Majorana neutrino masses, however, the cele- 



brated seesaw mechanism 
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29| assumes that lepton number is violated at a 



very large scale, probably at energies comparable to the grand unification scale. This "clas- 
sical" version of the seesaw mechanism, while automatically suppressing neutrino masses 
without the need for any small pre-factor, will unfortunately never be directly testable. 
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However indirect insight into the high-energy world might become possible, if weak scale 
supersymmetry is realized in nature. Indeed, starting from flavour diagonal soft super- 
symmetry (SUSY) breaking terms at some high-energy "unification" scale, flavour violation 
appears at lower energies due to the renormalization group running of the soft breaking 
parameters jsO]. If the (type-I) seesaw mechanism is responsible for the observed neutrino 
masses, the neutrino Yukawa couplings leave their imprint in the slepton mass matrices as 
shown first in 
potentially; 
3 



311 ] ■ Flavour off-diagonal entries in the neutrino Yukawas then can lead to 
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arge le pto n flavour violating lepton deca ys such as / 
I or /i — e conversion in nuclei 
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I]. In a similar spirit, if super- 
symmetry is discovered at a future accelerator such as the LHC, one can use measurements 
of masses and branching ratios of supersymmetric particles to obtain indirect information 



on the range of allowed seesaw parameters [40|, |41|, |42, |43|, |4J] . The most general supersym- 
metric seesaw mechanism has too many parameters to be predictive and thus can not be 
excluded by any measurements of lepton flavour violating (LFV) processes. Within the su- 
persymmetric version of the seesaw measurements of LFV observables outside the neutrino 
sector allow one to obtain valuable independent information about the seesaw parameters 



451]. There are two logical possibilities of how such LFV measurements might be useful, (a) 



Given the current incomplete knowledge on the light neutrino masses and angles, one could 
make some simplified assumptions about the right-handed neutrino sector. Then "predic- 
tions" for LFV observables as a function of the remaining unknowns for the left-handed 
light neutrinos result. Or (b) one could learn about the parameters of the right-handed 
neutrinos once the most important, but currently unknown light neutrino observables have 
been measured. While the second option might look more interesting, the time scale for 
making progress on m^, S13 or the Dirac CP phase 6 will be long. Worse still, the Majorana 
phases of the light neutrinos are unlikely to be ever reliably measured. Hence experimental 
information most likely will be incomplete and measurements of LFV observables will be 
useful to at least partially reconstruct the seesaw parameters. 

In this paper we study lepton flavour violating decays of the scalar tau as well as LFV lep- 
ton decays at low energies. We assume minimal supergravity (mSugra) boundary conditions 
and type-I seesaw as origin of neutrino masses and mixings. We compare the sensitivities of 
low-energy and accelerator measurements and study their dependence on the most important 
unknown parameters. LFV measurements at accelerators could be argued to be preferable 
to the low-energy LFV experiments for "reconstructing" seesaw parameters, since from a 
theoretical point of view they involve fewer assumptions. However, the absolute values of 
LFV stau decays and, for example, Br(/i ^ e + 7) depend very differently on the unknown 
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SUSY spectrum. Whether low-energy LFV or LFV at accelerators yields more insight into 
the seesaw mechanism can currently therefore not be predicted. 

While absolute values of LFV observables depend very strongly on the soft SUSY breaking 
parameters, it turns out that ratios of LFV branching ratios can be used to eliminate most 
of the dependence on the unknown spectrum. I.e., ratios such as, for example, Br(f2 — >■ 
^ + Xi)/Br(f2 A^ + Xi) cire constants (for fixed neutrino parameters) over large parts of the 
supersymmetric parameter space and therefore especially suitable to extract information 
about the seesaw parameters. We therefore study such ratios in detail, first in a useful 
analytical approximation and then within a full numerical calculation. 

The rest of this paper is organized as follows. In the next section, we will recall the 
basic features of the supersymmetric seesaw mechanism, mSugra and LFV in the slepton 
sector. Section UTTl then discusses analytical estimates for slepton mixing angles and the cor- 
responding LFV observables. In Sec. [IV] we present our numerical results before concluding 
in Sec. |Vl 



II. SETUP: MSUGRA WITH TYPE I SEESAW 



In order to fix the notation, let us briefly recall the main features of the seesaw mech- 
anism and mSugra. We will consider only the simplest version of the seesaw mechanism 
here. It consists in extending the particle content of the minimal supersymmetric standard 
model by three gauge singlet "right-handed" neutrino superfields. The leptonic part of the 
superpotential is thus given by 

W = Y^%HdE] + Yl%%N'r + M.iV^iV;. (2) 

where Ye and Yy denote the charged lepton and neutrino Yukawa couplings, while N'[ are 
the "right-handed" neutrino superfields with Mj Majorana mass terms of unspecified origin. 
Since the are singlets, one can always choose a basis in which the Majorana mass matrix 
of the right-handed neutrinos is diagonal Mpi. 

Note that LFV arises from supersymmetric as well as from gauge boson loop diagrams, 
for example slepton-gaugino exchange loops and W loops involving right-handed neutrino 
exchange. The former (SUSY-induced LFV) can be described by taking a basis where the Ye 
Yukawa coupling matrix is diagonal, its entries fixed by the observed charged lepton masses. 
This reduces the relevant physical parameters to a total of 21 . 

While in extended schemes like inverse seesaw 
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gauge-induced LFV is poten- 
tially sizeable, it is negligible in the simplest type-I seesaw model, due to the large values of 
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Mi required. Therefore, we focus on such intrinsically supersymmetric LFV, which can be 
well characterized Eq. ([2]) in the unbroken SU(2) limit. 

Different parametrizations for the simplest seesaw have been discussed in the literature. 
The most convenient choice for our calculation is to go to the basis where the charged lepton 
mass matrix is diagonal. We then have as parameters 9 mass eigenstates (3 charged leptons, 
the 3 light and the 3 heavy neutrinos). The remaining 12 parameters can be encoded in two 
matrices Vl and Vr, with 3 angles and 3 phases each, which diagonalize Y,^, 

% = V^Y^Vl. (3) 

The effective mass matrix of the left-handed neutrinos is given in the usual seesaw approxi- 
mation as 

2 

If one of the Mi eigenvalues of the matrix Mr goes to infinity (or the corresponding vector 
in Yy to zero) the corresponding eigenvalue of my {mi) goes to zero. Since the neutrino mass 
matrix is complex symmetric, Eq. (jlj) is diagonalized by 

tT 



= -^yJ ■ Mr' ■ (4) 



m,y = U ■ m^, ■ U. (5) 



Inverting the seesaw equation, Eq. (jll), allows one to express Fjy as [47 1 



Y, = V2—\JMr-R-^/^,-U^. (6) 

where rhy is the diagonal matrix with mj eigenvalues and R in general is a complex orthogonal 
matrix. Note, that in the special case R = \,Yy contains only "diagonal" products \/ Mimi. 
Note that in this approximation the 18 parameters in Y^ are reduced to 12, which are 
expressed as six light neutrino mixing angles and phases in the lepton mixing matrix t/, the 
3 light neutrino masses in m^ and the 3 heavy " right-handed" neutrino masses in \/~Mr. 

In the general MSSM, LFV off-diagonal entries in the slepton mass matrices are free 
parameters. In order to correlate LFV in the slepton sector with the LFV encoded in Yy 
one must assume some scheme for supersymmetry breaking. We will restrict ourselves here 
to the case of mSugra, characterized by four continuous and one discrete free parameter, 
usually denoted as 

mo, Ml/2, ^0, tan/5, Sgn(;u) (7) 

Here, m^ is the common scalar mass, M1/2 the gaugino mass and the common trilinear 
parameter, all defined at the grand unification scale, Mx — 2 ■ 10^^ GeV. The remaining two 
parameters are tan/5 = vjjlv p a nd the sign of the Higgs mixing parameter /x. For reviews 
on mSugra, see, for example [48|, |49|. 
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Calculable LFV entries appear in the slepton mass matrices, due to the nontrivial gen- 



eration structure of the neutrino Yukawa matrix in Eq. ([2]), as first pointed out in 31|. 



In order to determine their magnitude we solve the complete set of renormalization group 



equations, given in |33|, l50| • It is however useful for a qualitative understanding, to consider 



first the simple solutions to the renormalization group equations found in the leading log 



approximation [33|, given by 



(AM|)., = -^(3m^ + Al){Y}LY^),, (8) 



{AMl)i^ = 0, 

where only the parts proportional to the neutrino Yukawa couplings have been written. The 
factor L is defined as 

Equation ([HD shows that, within the type-I seesaw mechanism the right slepton parameters 
do not run in the leading-log approximation. Thus, LFV scalar decays should be restricted 
to the sector of left-sleptons in practice, apart from left-right mixing effects which could 
show up in the scalar tau sector. Also note that for the trilinear parameters running is 
suppressed by charged lepton masses. 

Note also that the LFV slepton mass-squareds involve a different combination of neutrino 
Yukawas and right-handed neutrino masses than the left-handed neutrino masses of Eq. (j4]). 
In fact, since (yJ-LFj,) is a hermitian matrix, it obviously contains only nine free param- 



eters 45|, the same number of unknowns as on the right-hand side of Eq. ([6]), given that 
in principle all 3 light neutrino masses, 3 mixing angles and 3 CP phases are potentially 
measurable ^. 

In an ideal world where all low energy paramaters, namely the 3 light neutrino masses, 
3 mixings and 3 CP violation parameters were known, the remaining parameters entering 
Eq. ([2]) could in principle be reconstructed by measuring all entries in (AM|)jj. This would 
determine the full set of 18+3 parameters which, to a good approximation, characterize LFV 
in the minimal type-I seesaw. In practice, however, there are two obstacles, (i) Calculability 
of (AM|)jj using Eq. assumes implicitly that there are no threshold effects near the 



^ In practice measuring the unknown angle 9i3 and the Dirac CP phase requires improved neutrino oscil- 
lation studies 5l| and will not be an easy task. Even if we are lucky to measure the overall neutrino 
mass scale in PPou experiments |5|], the Majorana phases contained in U are much harder to determine in 
practice. 
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unification scale which destroy the strict proportionality to the parameters mo and Aq 52|. 
In realistic grand unified theory models this might not be the case. And (ii) it is not realistic 
to assume that all entries in (AM|)jj can be measured with sufficient accuracy, since (a) 
the diagonal shifts (AM|)jj are very small compared to (M|)jj (nearly everywhere in the 
available parameter space) and (b) the determination of the phases require to measure CP- 
violating LFV observables. The latter does not seem to be a very realistic option either, 
since, as our numerical results show, one expects only rather low statistics to be available 
in measurements of LFV slepton decays. 



III. ANALYTICAL RESULTS FOR FLAVOUR VIOLATING PROCESSES 

In this section we present some general formulas describing lepton flavour violation within 
type-I seesaw schemes. We concentrate on the discussion of ratios of LFV branching ratios, 
since, as mentioned in the introduction, these are most easily connected to the seesaw 
parameters. As a first approximation we adopt the mass insertion approximation, neglecting 
left-right mixing in the slepton mass matrix and taking the leading-logs (see below). We will 
demonstrate the reliability of our analytical estimates in the next section, where we perform 
a full numerical calculation of the various LFV branching ratios, which does not rely on any 
of the approximations discussed in this section. 



A. General formulas 



The charged slepton mass matrix is a (6,6) matrix, containing left and right sleptons. 
Here we concentrate exclusively on the left-slepton sector. Taking into account the discussion 
given in Sec. [Til this is a reasonable first approximation, as can be seen from Eq. ([8]). The 
left-slepton mass matrix is diagonalized by a matrix R\ which in general can be written 
as a product of three Euler rotations. However, if the mixing between the different flavour 
eigenstates is sufficiently small, i?' can be approximated as 

/ 1 Opr, 6z~ \ 



V 



-9, 



9 



(10) 



1 J 



an approximation that corresponds to that employed in the mass-insertion method SlJ. In 
this small-angle approximation each angle can be estimated by the following simple formula 



9- ■ ~ 



(AM2 



(AM|),,-(AM|),,' 



(11) 
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LFV decays are directly proportional to the squares of these mixing angles, for example 
Br(/x — s> e + 7) ~ {GefiY if all angles are small. 

Within mSugra ratios of LFV branching ratios can then be used to minimize the depen- 
dence of observables on SUSY parameters. Consider the case of LFV decays which involve 
only one generation of sleptons, for example Br (7-2 e + Xi) ciiid Br(f"2 ^ yU + x?)- Taking 
the ratio of these two decays 



i?r(f2^e + x?) /^ef\2 /(AM? 



5r(f2^/i + X?) V^/if^ V(AM? 



J 13 



j23 



(12) 



i.e. one expects that (a) all the unknown SUSY mass parameters and (b) the denomina- 
tors of Eq. f|TT]) cancel approximately. The latter should happen practically everywhere in 



mSugra parameter space since (M?)ee — (M?)^^. This straightforward observation forms 
the basis for our claim that ratios of branching ratios are the theoretically cleanest way to 
learn about the unknown seesaw parameters. Numerically we have found, that relations 
similar to Eq. 0121) hold also for ratios of observables involving decaying particles of different 
generations, such as the low-energy ratio Br(/i — > e + 7)/Br(r — e + 7). 

To calculate estimates for the different ratios of branching ratios we therefore define 



I (AM?) 



kl\ 



where the observable quantity is (?"fc;)^. Of course, only two of the three possible combinations 
that can be formed are independent. For example, Br(/i — e+7) /Br(r — > e-l-7) ~ {rlD'^xTZ. 
Here, 7^ is a correction factor taking into account the different total widths of the muon and 
the tau, 7^ = F^/F^ ^. 

In the leading-log approximation the off-diagonal elements of the charged slepton mass 
matrix are proportional to (AM|)jj oc ((y^)^L(y ))^^.. Using the parametrization for the 
Yukawa couplings of Eq. the entries in (AM|)jj can be expressed as 

(AM^),,- (X U,^U*^y/^^Rl^RkpMk log (^^^ . (14) 

We can now rewrite Eq. ( fT4|) in terms of observables which are more directly related to 
experiments. In the standard parametrization for the leptonic mixing matrix U is completely 



The inclusion of this factor (and similar corrections for the other low-energy LFV decays) is necessary, 

4) 



since (7"^;)^ relate really partial widths, whereas the measured quantity is usually the branching ratio. 
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analogous to the CKM matrix and can be written as 

12C13 ^136 ^"^^ 

(15) 



^ C12C13 S12C13 Si3e 



-S12C23 - Ci2S23Sl3e*'' C12C23 - Si2S23Si3e*'^ S23C13 
y S12S23 — Ci2C23Si3e*'' — C12S23 — Si2C23Si3e*'' C23C13 J 



where we assumed strict unitarity and neglected the Majorana phases 8|, because they do 
not affect lepton number conserving processes such as the LFV decays we are concerned 
with here. 

Given that neutrino oscillation experiments fix two mass-squared splittings, we can re- 
express the three light neutrino masses in terms of one overall neutrino mass scale and the 
measured quantities AttIq and Am^, where Am^ (Am^) is the solar (atmospheric) mass- 
squared splitting. We will refer to the case of m\ = {rris = 0) as strict normal (inverse) 
hierarchy. This choice has the advantage that in both cases S12 = sin^Q and S23 = sin6'A. 
Equation (fT4l) can then be written in terms of the measured neutrino angles su and S23, 
the measured neutrino mass-squared splittings, plus the so far unknown overall neutrino 
mass scale rrii, and the reactor neutrino angle S13 = sr. If the latter were measured, one 
could extract information on the right-handed neutrino mass scale and/or the matrix R from 
Eq. iHM . Conversely, we could learn about m,^ and S13 from measurements of LFV decays, 
making some assumptions about the scale Mr and the possible textures of the Yukawa 
couplings that determine M/j and R. 



B. Degenerate right-handed neutrinos 

In this subsection we will assume that the three right-handed neutrinos are degenerate. 
This simplifying ansatz allows us to study the sensitivity with a single mass-scale parameter 
associated with the neutrino mass generation via type-I seesaw mechanism. This ansatz 
can be theoretically motivated in the framework of some flavour symmetries, for example 
A4 53|. In the special case that the matrix R is real, Eq. (fT4|) reduces to 

(AM|)^2 ^ C12C13 (-S12C23 - Ci2S23Si3e"*'^) 2l (16) 

+ ■5i2Ci3 (C12C23 — Si2S23Si3e *'') Z2 + 5230135136 

(AM|)^3 OC C12C13 (S12S23 - Ci2C23Si3e"*'') Zl 

+ 512013 ( — C12S23 — Si2C23Si3e *'') Z2 + £230135136 Z3, 

(AM|)23 OC (S12S23 - Ci2C23Si3e^*'^) (-S12C23 - Ci2S23Si3e*'') Zl 

+ (-C12S23 - Si2C23Sl3e"*'') (C12C23 - Si2S23Si3e*'') Z2 

+ S23C23CI3Z3, 
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where 



rriiMi log 



M 



X 



Mi 



(17) 



For this degenerate right-handed neutrino ansatz the combination Mjlog(^^) becomes an 
overall factor, which can be taken out from Eq. ffUl) . since it cancels upon taking ratios. 
I.e. for degenerate right-handed neutrinos one may simply make the replacement Zi —>■ rrii 
in Eq. (fT6l). 

As a starting approximation for the following estimates, let us assume that the lepton 
mixing matrix has the exact tribimaximal (TBM) form 54j ] 



U = U^ 



TBM 







V 3 

1 1 1 
V3 V2 



(18) 



As is well-known, Eq. ( 1181) is an excellent first-order approximation to the measured neutrino 
mixing angles 4] • With this assumption the ratios of the off-diagonal elements of the charged 
slepton mass matrix are simply given by 



,12 
13 



.12 
23 



(19) 



.13 
23 



2(m2 — nil] 



|3m3 — 2m2 — mi\ 



As Eq. ( p^Qj) shows, and depend on mass squared splittings and on the overall neutrino 
mass scale, i.e. also on the unknown neutrino mass hierarchy. In the case of strict normal 
hierarchy (SNH, mi = 0) 



,12 



13 



^"23 ~ ^"23 



2v^ 



(20) 



where a = ^'"^ 



I Ami I 



while for the case of strict inverse hierarchy (SIH, m-^ = 0) 



,12 



13 



^"23 ~ ^"23 



2(1 -vO^) 



2 + 

Finally, for quasi-degenerate (QD) neutrinos, defined as a/ Am\ <^ m,^, one finds 



(21) 



/>,12 /V.13 

'23 ~ ' 



2a 
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where cta is the sign of the atmospheric mass splitting 
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FIG. 1: Ratio versus the neutrino mass ?ni in eV. The light/yellow (dark/violet) band is for 

the case of normal (inverse) hierarchy. The width of the band indicates the uncertainty due to the 
currently allowed 3 a C.L. ranges for Am^ and AttIq. The calculation assumes exact tribimaximal 
mixing for the left-handed neutrinos. 



Note that cr^ equals +1 (—1) for normal (inverse) hierarchy. Thus QD neutrinos with 
normal (QDNH) or inverse hierarchy (QDIH) lead formally to different results. However, 
this difference is numerically not relevant, once uncertainties are taken into account. 

Figure [T] shows the ratio (^gf )^ versus the neutrino mass mi in eV for normal (inverse) 
hierarchy. The figure demonstrates the importance of the absolute neutrino mass scale for 
(^23)^- ■'■^ most general case one must use Eqs. f|T3|) and f|T6l) . However, for S13 = 
the explicit dependence of the ratios of the off-diagonal elements of the charged slepton 
mass-squared matrix on the other neutrino angles matrix follow rather simple expressions 



,12 

13 



„12 
23 



„13 
23 



m2 — mi 



C23 

■§23 
1 

■Sl2Ci2-j 

S23 |"^3 

1 

SI2C12-, o o 

C23 I"i3 - cfama - sl^mi 



(24) 



■ C?2"^2 

m2 — mi 



sj^mil' 



Figures [2] and [3] show the dependence of the square ratios {r]^iY as a function of 5^3 for the 
different extreme cases of SNH and SIH as well as QDNH and QDIH, for two choices of the 
Dirac phase 5 = 0, vr. These ratios (r^])^ depend strongly on sl^. Note from Eq. f|T6|) that for 
tan^ 6a = 1, (^"23)^ (^"23)^ invariant under exchange of 5 = ^ 5 = tt. If tan^ Oa 7^ 1, 
this symmetry is broken, but always one of the two ratios rgf and is guaranteed to be 
non- vanishing regardless of the value of S13. A non-zero measurement of both ratios would 
therefore in principle contain information on both S13 and 6 (if right-handed neutrinos are 
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degenerate). 




, „-4 
10 



43 



^3 



FIG. 2: Square ratios {r\^)^ (blue line, dotted line), (^23)^ (red line, dashed line) and (rgf)^ (g 
line, full line) versus 5^3 for SNH (upper panels), SIH (lower panels) for 5 = (left panels) 
6 = IT (right panels). The plots assume that the heavy neutrinos are degenerate. The other 
neutrino parameters have been fixed to their b.f.p. values. Note from Eq. (jl6p . that for tan^ 9a 



[Teen 
and 

light 
= 1, 



{rliY and (r. 



23)^ are symmetric under the exchange of 5 = ^ 5 



TT. 



Figure [3] shows also that the cases QDNH and QDIH are also symmetric under the 
simultaneous exchange of(5 = 0^5 = 7r and QDNH ^ QDIH, for the case of tan^ = 1- 
This symmetry is broken in all cases for tan^ 7^ 1, as seen from the numerical values given 
in Tables |T] and [Tll Tables [T] and |TT] show numerical values for r]^i for the extreme cases of 
SNH, SIH, QDNH and QDIH for various different choices of neutrino parameters. In the 
rows labeled as TBM, we have used the TBM values for 612 and ^23and the neutrino mass 
splittings have been fixed to their best-fit values taken from Ref. In the rows labeled 
as 3(7, we take into account the experimentally allowed 3a ranges for neutrino oscillation 
parameters: sj^ = 0.26 - 0.40, 5^3 = 0.34 - 0.67, Am^ = (7.1 - 8.3) x 10-^^ ^^^d 
Am\ = (2.0 — 2.8) x 10~^eV^. In the first column, Oi^ has been fixed to its TBM value 
(■§13 = 0), while in the second and third columns S13 has been fixed to s^^^, which is the 
experimentally allowed maximum value: (s^^^^)^ = 0.050 at 3cr C.L. In the second column, 
the Dirac phase is fixed to 5 = 0, while in the third column 6 = n. Note that, as already 
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FIG. 3: As Fig. [2 but for the limit of quasi-degenerate light neutrinos. Square ratios (rjl) (blue 



12\2 



(red line, dashed line) and (green line, full line) versus s^g for QDNH 



line, dotted line), {r^j, 
(upper panels), QDIH (lower panels) for 5 = (left panels) and 5 = vr (right panels). 



mentioned, these estimates are valid in the small mixing limit and hence these values are 
indicative only. 



C. Right-handed neutrinos strongly hierarchical 

One can consider the case of degenerate right-handed neutrinos to be just one extreme 
limit in a continuum of possibilities. The opposite extreme case would than be to assume 
right-handed neutrinos are strongly hierarchical. Note that here we make the important 
assumption that the matrix R is the identity. 
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TABLE I: The parameters r^^^ are given for several values of neutrino oscillation parameters. SNH 
and SIH are strict normal and strict inverted hierarchy of neutrino masses, respectively. Rows 
labeled as TBM assume the TBM values for 9^2 and (?23 and the neutrino mass splittings have 
been fixed to their b.f.p. values taken from [J]. Rows labeled as 3(T take into account current 
allowed 3cr ranges of neutrino oscillation parameters. In the first column, ^13 has been fixed to its 
TBM value (si3 = 0), while in the second and third columns S13 has been fixed to its maximum 
allowed value: [sff^Y = 0.050 at 3(7 C.L. and the Dirac phase is fixed to (5 = and 5 = vr, 
respectively. 



1. Dominant Mi 

If Ml is the heaviest mass eigenvalue, the leading terms for the off-diagonal slepton masses 
are (in case mi 7^ 0) 

(AM|)^2 ^ Ci3Ci2(Si2C23 + Si3e-^^'^Ci2S23) (25) 
(AM|)^3 OC Ci3Ci2(Si2S23 - Si3e-*''ci2C23) 

(AM|)23 OC 312523023 - Si3Si2Ci2(e"*''c23 - e*''s23) - 523012523023 

For the special case of 513 = 0, the ratios simplify to rj? = — , r,? = — ^nd r\l = . 
Note the large difference in the numerical values compared to the case of degenerate right- 
handed neutrinos. Here, for example for S13 = one finds (^21)2 = 4, whereas in the case of 
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TABLE II: The parameters r^-'^ are given for several values of neutrino oscillation parameters. QD 
stands for the quasi-degenerate limit, while NH (IH) indicate that the neutrino hierarchy is normal 
(inverse) . The neutrino parameters have been varied in the same way as in Table HI 



degenerate right-handed neutrinos one obtains (^21)2 = 0.017 [best fit point (b.f.p.) values 
for AttIq and Am\]. For nonzero values of sis Fig. H] shows that [r^j^iY depend to a much 
lesser degree on S13 than for the case of degenerate right-handed neutrinos. Especially, note 
that for the case of Mi dominance considered here none of the {r]^i)'^ vanish in the allowed 
range of S13. Numerical values for extreme values of S13 are summarized in Table UTTl 



2. Dominant M2 

If M2 is the heaviest mass eigenvalue, the dominant terms for the off-diagonal slepton 
masses are 

(AM|)^2 °^ Ci3Si2(ci2C23 - Si3e-^^Si2S23) (26) 
(AM|)^g OC Ci3Si2(ci2S23 + Si3e"*'^Si2C23) 

(AM|)23 OC C22S23C23 + ■Si3Si2Ci2(e"*''c23 - e^'^Saa) - S23S22S23C23 
For the special case of si^ = 0, the ratios simplify to rl? = — , r,? = and ri? = . 

f ±01 f J Li S23^ ^•^ C12S23 ■^•^ C12C23 

Here, for example, for S13 = one finds (^21)2 = 1, whereas for the case of Mi being 



15 




FIG. 4: Square ratios (rjl)^ (blue line, dotted line), (^23)^ (red line, dashed line) and (^23)^ (green 
line, full line) versus S13 for 5 = (left panel) and 5 = n (right panel) for the case of Mi being 
dominant. The remaining neutrino parameters have been fixed to their b.f.p. values. 

dominant this quantity is expected to be (^"23)^ = 4. Figure [5] shows the (r^j)^ as function 
of s^3 for the M2 dominance case. Again the dependence on S13 is weaker than in the case 
of degenerate right-handed neutrinos. As in the previous case (r^j)^ never vanishes in the 
allowed range of sl^. Finally, the numerical values also differ from the ones found for the 
case of Ml dominance. A summary of numerical values for extreme values of S13 is given in 
Table Iml 




FIG. 5: Square ratios (r||)^ (blue line, dotted line), (^2!)^ (red line, dashed line) and (^23)^ (green 
line, full line) versus S13 for 6 = (left panel) and 6 = it (right panel) in the case where M2 is 
dominant. The other neutrino parameters have been fixed to their b.f.p. values. 
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3. Dominant M3 

If terms proportional to M3 give the leading contribution one finds 

(AM?) ^2 ^ si3e"''ci3S23 (27) 

(AM?) ^3 (X Si3e-'^Ci3C23 
(^^1)23 °^ ^?3^23C23 

For the special case of S13 = 0, one finds that = r^^ = 0, otherwise both ratios are 
proportional to S13. These numerical values allow us to distinguish the M3 dominance case 
from the previous hierarchical cases already discussed. Numerical values for extreme values 
of Si3 are summarized in Table lllli 



IV. NUMERICAL RESULTS 

The analytical results presented above allow us to estimate ratios of branching ratios for 
LFV decays. For absolute values of the branching ratios, as well as for cross-checking the 
reliability of the analytical estimates, one must resort to a numerical calculation. In this 
section we present results of such a numerical calculation. All results presented below have 
been obtained with the lepton flavour violating version of the program package SPHENO [55 
For definiteness we will present results only for the mSugra "standard points" SPSS [5 



and SPSla' 



57| . taken as reference examples. However, we have checked with a number 



of other points that our results for ratios of branching ratios are generally valid. SPSla' 



57( 1 is a typical point in the "bulk" region for SUSY dark matter. It is a slightly modified 



version of the original SPSla point of [56|], which gives better agreement with the latest 
constraints from cold dark matter abundance. It has a relatively light slepton spectrum, i.e. 
left sleptons around 200 GeV. SPSS 56| is a point in the co-annihilation region for SUSY 
dark matter. Left sleptons in this point are heavier than in SPSla', i.e. have masses around 
S50 GeV. We have chosen these two points to show the complementarity between low-energy 
searches for LFV and LFV scalar tau decays at the LHC, see also the discussion below. 

Our numerical procedure to fit the neutrino masses is as follows. Inverting the seesaw 
equation, see Eq. (jlj), one can get a first guess of the Yukawa couplings for any fixed values of 
the light neutrino masses and mixing angles as a function of the corresponding right-handed 
neutrino masses. We then run numerically the renormalization group equations taking into 
account all flavour structures in matrix form. We integrate out every right-handed neutrino 
and its superpartner at the scale corresponding to its mass and calculated the corresponding 
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TABLE III: The parameters r^^ are given for several values of neutrino oscillation parameters. 
Each row labeled as Mj is calculated assuming the contribution from neutrino with mass Mj is 
dominant. Neutrino oscillation parameters have been varied as in Table HI Notice that the row for 
dominant M3 gives the same numerical result for the Dirac phase 5 = and 6 = tt. 



contribution to the dimension-5 operator which is evaluated to the electroweak scale. This 
way we obtain the exact neutrino masses and mixing angles for this first guess. The difference 
between the results obtained numerically and the input numbers is then minimized in an 
iterative procedure until convergence is achieved. As is well-known neutrino masses and 
mixing angles run very little if physical light neutrino masses are hierarchical 58|. Thus, 
barring the exceptional case where neutrinos become very degenerate, one usually reaches 
numerical convergence very fast. For degenerate left neutrinos convergence from first guess to 
exact results can be slow, especially for relatively large values for the right-handed neutrino 
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masses, which require larger Yukawa couphng constants. In this case we used a numerical 
fit procedure 59| based on the program MINUIT 

In the following two subsections we present numerical results first for the case of degener- 
ate right-handed neutrinos, then for the case(s) of very hierarchical right-handed neutrinos. 
We have checked numerically that, as expected from Eq. ([8]), right-sleptons have small 
branching ratios for LFV final states. Thus, the discussion concentrates on the decays of 
the "left" staus f2 ^ tl. 



A. Degenerate right-handed neutrinos 

In this subsection we still adopt the simplifying ansatz that R = 1, see Eq. ([6]). Two 
examples for hierarchical light neutrinos are shown in Fig. P and Fig. [71 Figure E] has the 



mSugra parameters fixed to the standard values SPSla' |56|, 157|, while Fig. [7] corresponds to 
SPSS 561]. The neutrino oscillation data are fitted for the strict normal hierarchy (SNH) case 
where mi = with exact tribimaximal mixing. The plot on the left panel shows low-energy 
lepton fiavour violating decay branching ratios for /j ^ /j + 7 and /j — 31 j, while the one on 
the right panel gives LFV stau (f2) decay branching ratios as a function of the right-handed 
neutrino mass scale Mi = Mr. 

As expected, all LFV processes show a strong dependence on M/j. This can be straight- 
forwardly understood from Eqs. (jlj) and ([8]). Keeping the light neutrino masses constant 
AM| are proportional to M^ log Mr, thus all LFV branching ratios grow as (MRlogM^)^. 
As the figures show, as long as Mr is not too large, all lepton fiavour violating processes 
show the same dependence on Mr. Ratios of branching ratios follow very nicely the corre- 
sponding analytically calculated ratios for {r]^i)^, once the corresponding correction factors 



are taken into account for the low-energy observables. As is well known 35|, l39|, |60|, for 
most parts of the mSugra parameter space one expects 

,2 



Br(/i^/, + 7) SvrV ^^mf/ 4 A 



(28) 



thus the photonic penguin diagram dominates the three- lepton decay modes 31 j. 

Figures E] and [7| do indeed confirm the validity of this approximation. Only at large 
values of Mr one observes some deviations from the analytical estimates. The reason for 
this departure is that in this parameter range the small-angle approximation no longer 



^ Minimization package from the CERN Program Library. Documentation can be found at 
http : / / cernlib . web . cern . cli/cernlib/ 
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holds, as can be seen from the absolute values for the decay Br(f2 ^ fi + Xi)) which can 
reach more than 10 % for M/j > 10^^ GeV. However, Figs. E] and [7] also show how the 
LFV f2 decays are strongly constrained by low energy data. For the degenerate right- 
handed neutrino case shown here (and for S13 = 0), independent of the mSugra parameters 
Br(/x — s> e + 7) is the most important constraint. Applying the current experimental limit on 
Br(/x — s> e + 7) of Br(/i — > e + 7) < 1.2 x 10^"*^^ {6l|, the branching ratio for Br(f"2 — > /i + Xi) 
is expected to lie below 10"'^ for SPSla', whereas it can reach several percent in case of 
SPSS. Note that in the range of M^j not excluded by the limit on Br(/i ^ e + 7) the ratio 
Br(f2 ^ e + Xi)/Br(f2 — ^ // + Xi) follows very well the analytical estimate of Eq. ( l20l) . 
The huge difference in the upper limit for Br(f2 fi + Xi) when going from SPSla' to 
SPSS can be understood from the fact that both left-sleptons as well as (lightest) neutralino 
and chargino are approximately a factor of two heavier for SPSS than for SPSla'. Since 
Br(/i e + 7) oc ^/mgjjgy |23| one expects Br(/i — > e + 7) to be a factor of more than 
several hundred lower for SPSS than for SPSla'. 




FIG. 6: Branching ratios for li — > + 7 and li 
(right) for the standard point SPSla' versus M/j, assuming degenerate right-handed neutrinos. 
Neutrino oscillation parameters have been fixed to the best fit values for Amg and Ainj^, with 
exact tribimaximal neutrino angles. We also set mi = 0. The coloured region in the right-side plot 
is excluded from the current experimental limit on Br(/x ^ e + 7). Thus, one expects for SPSla' 
only very small branching ratios for LFV scalar tau decays. (Compare to Fig. ([7]). 

The strong dependence of Br(/i e + 7) on the supersymmetric mass spectrum is also 
seen in Fig. |H1 where we plot Br(/i — > e + 7), Br('f2 — /i + Xi) Br(f2 e + Xi) versus 
the mass of f2, for light neutrino parameters as before and a fixed value of Mr = S x 10^^ 



31 j (left) and f2^e + Xi and f2^n + Xi 



20 




10'° 10" 10" 10'° 10'° 10" 10'^ 10" 10" 10'= 10'= 

M, (GeV) M, (GeV) 



FIG. 7: Same as Fig. [6l but for the mSugra standard point SPS3. In this point the constraints on 
the LFV f2 decays from the upper hmit on /_f ^ e + 7 are much less severe than for SPSla'. As a 
result Br(f2 — > fi + Xi) could be as large as several percent with all low-energy constraints fulfilled. 
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FIG. 8: Branching ratios as function of scalar tau mass. The full line (red) is 10^^- Br(;U — > e + 7), 
the dashed line (blue) Br(f2 ^ /U + Xi) and the dot-dashed line (green) is Br(f2 — > e + Xi)- Data 
calculated for SPSla with parameters varied along the "slope". Note that SPSla is used in this 
plot instead of SPSla', since for SPSal' no slope is given in {s?!. Right-handed neutrino mass is 
fixed to Mji = 3 x 10^'^ GeV. The black line is the current upper limit on Br(/i — > 67). While 
SPSla with Mji = 3 x 10^^ GeV is excluded by Br(/_f 67), for slightly heavier slepton masses the 
low-energy constrained can be evaded, having at the same time sizeable lepton flavour violating 
slepton decay branching ratios. 
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GeV. Here, the parameters for the point SPSla have been varied around the slope given 
in Ref. 56|]. Note that Br(/i ^ e + 7) drops below the current experimental limit for ruf^ 



larger than about 250 GeV. In contrast, the f2 LFV decay branching ratios increase for 
increasing rrif^. This is due to the fact that left-sleptons become more degenerate when 
mo is increased along the slope for SPSla. The more degenerate sleptons are, the larger 
the resulting LFV parameters, for given light neutrino parameters. Note, however, that the 
ratio Br(f2 — > e + Xi)/ Br(f2 —* fi + Xi) remains constant in agreement with the analytical 
estimate, as long as Br(f2 ^ fi + Xi) is smaller than a few percent. Again this reflects the 
fact that the small-angle approximation is valid only for small branching ratios in the LFV 
decays. 

We have also checked numerically the reliability of our analytical calculation for the case 
of mi ^ 0. An example is shown in Fig. [91 Here we have fixed the mSugra parameters to 
the standard point SPSla', the right-handed neutrino mass scale to Mr = 5 X 10^^ Qgy, the 
light neutrino mixing angles to the TBM values, Am^ and Amg to their b.f.p. values and 
we have calculated Br(f2 — * e + Xi)/Br(f2 ^ fi + Xi) ^s a function of the lightest neutrino 
mass. As shown in Fig. [9] the value of this ratio obtained within a full numerical calculation 
follows very closely the central value given in Fig. [H as expected (here we assumed the case 
of normal hierarchy). 




FIG. 9: Numerically calculated Br(f2 ^ e + Xi)/Br(f2 ^ fi + Xi) for the standard point SPSla' 
versus lightest neutrino mass for the case of normal hierarchy, (compare to Fig. [1]). 
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B. Hierarchical right-handed neutrinos 

Now we turn to the extreme case of very hierarchical right-handed neutrinos. Again 
our goal is to check the reliability of the analytical calculation for this case. In all figures 
presented in this subsection we have taken two of the three right-handed neutrino masses to 
be constant at Mf> = 10^° GeV and varied the remaining third right-handed neutrino mass 
in the ranges given in the figures. In all cases we have fixed the neutrino angles to the TBM 
values, Amj^ and AniQ to their best-fit values and assumed normal hierarchical neutrinos. 
The remaining free parameter nii is given in each figure. 




10'° 10" 10" 10" 10'° 10" 10" 10' 

M, (GeV) M, (GeV) 



FIG. 10: Branching ratios for + 7 ™d li 31 j, as a function of Mi for constant M2 = 

M3 = 10^*^ GeV and nii = 0.001 eV (left) and for mi = 0.1 eV (right). mSugra parameters have 
been fixed to SPSla'. 

Figure [TO] shows LFV lepton decays as a function of Mi for mi = 0.001 eV (left) and for 
mi = 0.1 eV (right) for the mSugra parameters fixed at SPSla'. For mi = 0.001 eV, the 
curves are not monotonous functions of Mi. In fact, in the left figure only for Mi > 10^^ 
GeV do the different branching fractions follow the analytical estimates of Eq. ( 125|) . This is 
due to the fact that the different contributions of the Mi to AM|_ scale like mjMj log Mj, 
i.e. Ml becomes dominant in the expressions for the AM|_ only if Mi/Mj ^ mj/mi. This 
is confirmed by the figure in the right panel, for which mi = 0.1 eV has been chosen. Here, 
the contribution from Mi to the AM? is indeed the dominant one for Mi > (few) xlO^° 
GeV. 

Figure [11] shows branching ratios for f2 e(/i)-|-Xi as a function of Mi for the two mSugra 
points SPSla' (left) and SPSS (right). Again the region excluded by the current upper limit 
on Br(/i ^ e -|- 7) is indicated. Ratios of the LVF slepton decays follow the analytical 
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FIG. 11: Branching ratios for f2 e(/z) + Xi SiS a function of Mi for constant M2 = M3 = lO^'^ 
GeV for SPSla' (left) and SPS3 (right). 

estimate very well everywhere in the region allowed by the upper limit on Br(/i e + 7). 
One observes, as is the case also for degenerate right-handed neutrinos, that for SPSla' the 
absolute values for the LFV branching ratios are too small to be observable, whereas for the 
mSugra point SPSS much larger values for LFV scalar tau decays are allowed. Note that 
Br(f2 — > e + Xi) is larger than Br(f2 — /i + Xi) for ^1 dominance, in contrast with the case 
of degenerate right-handed neutrinos. 
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FIG. 12: Branching ratios for li Ij + ^ and li — > 2>lj (left) and LFV stau decays (right), for the 
standard point SPSS as a function of M2 for constant Mi = M3 = 10"'^'^ GeV. 

Figure [12] shows branching ratios for li Ij + 'y and /j 31 j (left) and LFV stau decays 
(right), for the standard point SPSS as a function of M2. As in Fig. [101 the left panel 
illustrates that only for M2 > 10^^ GeV the contribution from M2 to the LFV mixing angles 
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is dominant. For M2 > 10^^ GeV the ratios of branching ratios follow the expectation of 
Eq. (1261) . LFV scalar tau decays as large as 1 % are allowed in this example. Note also that 
Br(f2 — e + Xi) = Br(f2 ^ fi + Xi) for M2 dominance and TBM neutrino angles. 
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FIG. 13: Branching ratios for li Ij + J and /j — > 31 j (left) and LFV stau decays (right), for the 
standard point SPSS as a function of M3 for constant Mi = M2 = W^^ GeV. 

Finally, Fig. [13] shows branching ratios for li ^ Ij + '~f and U — > 3lj (left) and LFV stau 
decays (right), for the standard point SPSS as a function of M3 fixing S13 = exactly. This 
implies that all final LFV states involving electrons are tiny, as is expected from Eq. fl27j) . 
Therefore for S13 = and M3 dominance there is no constraint from the upper limit for 
Br(/x — i> e + 7). Once S13 is nonzero branching ratios for LFV final states involving electrons 
also become nonzero and proportional to sf^. 

In summary this section demonstrates that the full numerical calculation confirms the 
analytical estimates presented above. Absolute values of the LFV branching ratios for lepton 
decays are sensitive functions of the unknown SUSY spectrum. For light sleptons, usually 
the constraint from the non-observation of Br(/i e + 7) makes the observation of LFV 
stau decays more likely when M3 gives the leading contribution to the LFV slepton mixing 
angles and S13 is close to zero. In this case LFV stau branching ratios may exceed 10%, as 
seen in Fig. [131 LFV stau branching ratios exceeding a percent are also possible for SPS3 
for hierarchical right-handed neutrinos and Mi and M2 dominance, as seen in Figs. [11] and 
W2[ but not for the SPSla' case. Similarly, for the case of degenerate neutrinos, LFV stau 
branching ratios can exceed a few percent, as seen in Figs. U\ especially for heavier sleptons, 
say 250 — 300 GeV, where the Br(yU — > e + 7) is smaller than the experimental limit and 
hence does not place a restriction, as seen in Fig. [8] 

Finally we note that we have expressed our results in terms of branching ratios. To get 
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a rough idea on the observabihty of the signal, one has also to consider cross sections and 
backgrounds. For the signal itself one would have to work out a detailed set of cuts to 
suppress background which is clearly beyond the scope of the present work. However, after 
applying basic cuts used for SUSY signals 62| one can estimate the cross sections for f2 



production. Using PYTHIA 6.4 63|] we find for the sum of all (Drell-Yan) cross sections 126 
fb (25 fb) and 31 fb (3 fb) for fo in cascade decays in the case of SPSla' (SPS3). Based 
on Monte Carlo analysis 6^, l65|] it has been shown that lepton flavour violation can be 
observed in dilepton invariant mass spectra within SUSY cascade decays. There the largest 
SM background is due to ti production. There is also SUSY background due to uncorrelated 
leptons stemming from different squark and gluino decay chains. The di-lepton spectra can 
provide a distinct signal of lepton flavour violation, namely the appearance of double peaks 



66l | due to the fact that not only one but two or more sleptons can contribute to these 



spectra. In case of Drell-Yan processes the main background will be W production. To 
show more clearly the observability of such LFV signals a detailed Monte Carlo study would 
be necessary. This, however, is beyond the scope of the present paper. 

We have shown results only for two "standard" mSugra points. However, as mentioned 
above, we have checked with a number of other points that ratios of branching ratios to a 
good approximation do not depend on the mSugra parameters. For absolute values of the 
branching ratios in general a heavier slepton spectrum leads to smaller LFV rates at low 
energy and larger LFV branching ratios at the LHC become possible, see also Fig. [SI Heavier 
sleptons, on the other hand, will lead to lower Drell-Yan production cross section, such that 
stau production will be dominated by cascade decays, the exact number of events depending 
on the details of the SUSY spectrum. We plan to do a more detailed, quantitative study of 
absolute event rates over all of mSugra space in the future. 



V. CONCLUSIONS AND OUTLOOK 

We have calculated lepton flavour violating processes both in LFV decays of the fi and the 
r leptons, as well as branching ratios for LFV stau decays in the supersymmetric version of 
the minimal type-I seesaw mechanism with mSugra boundary conditions. We have limited 
ourselves to the study of a few standard mSugra points, ratios of LFV branching ratios are 
independent of this choice and therefore an interesting instrument to study the unknown 
seesaw parameters. 

We have shown that the LFV branching ratios for lepton decays are sensitive functions 
of the unknown SUSY spectrum. For light sleptons, the non-observation of Br(/i — > e -|- 7) 
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places an important constraint on the observability of LFV stau decays. The most favorable 
case is when right-handed neutrinos are hierarchical, with M3 giving the leading contribution 
to the LFV slepton mixing angles and S13 close to zero. In this case LFV stau branching 
ratios may exceed 10% or so, see Fig. [T31 LFV stau branching ratios exceeding the percent 
level may also occur for hierarchical right-handed neutrinos with Mi or M2 dominance for 
the SPSS reference point, but not for the SPSla' case, see Figs. [TT] and [T2l Similarly, for the 
case of degenerate neutrinos, LFV stau branching ratios can exceed a few percent, as seen 
in Figs. U\ especially for sleptons heavier than 250 GeV or so, as seen in Fig. [HI 

Notice that the above results rely crucially on an important simplifying assumption about 
the right-handed neutrino spectrum. For example, for degenerate right-handed neutrinos 
they require that R be real, while for hierarchical right-handed neutrinos they hold when 
R = 1. This simplification allows one to calculate LFV decays of leptons and of the scalar 
tau as a function of low-energy neutrino parameters. However the use of this assumption 
should be critically scrutinized. We plan to come back to this issue in a future publication. 
Once an improved experimental determination of mi and S13 become available from future 
double beta decay and neutrino oscillation studies at reactor and accelerators, one could 
start "learning" about the right-handed neutrino sector, once the correct SUSY breaking 
scheme has been identified and provided that the SUSY breaking scale is above the lepton 
number breaking scale. 
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